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ABSTRACT 

Context. Studies of the infrared (IR) emission of cosmic sources have been proven essential to constrain the evolutionary history 
of cosmic star formation and that of nuclear black hole gravitational accretion, as the bulk of such events happens inside heavily 
dust-extinguished media. 24 pm). 

Aims. The Spitzer Space Telescope has recently contributed a large set of data to constrain the nature and cosmological evolution 
of infrared source populations. We exploit in the present paper a large homogeneous dataset to derive a self-consistent picture of IR 
emission based on the time-dependent A e ff=24, 15, 12 amd 8 pm monochromatic and bolometric IR luminosity functions (LF) over 
the full < z < 2.5 redshift range. 

Methods. Our present analysis is based on the combination of data from deep Spitzer surveys in the VIMOS VLT Deep Survey 
(VVDS-SWIRE) and GOODS areas. To our limiting flux of 5 24 = 400 pJy our derived sample in VVDS-SWIRE includes 1494 
sources, and 666 and 904 sources brighter than 5 24 = 80 pJy are catalogued in GOODS-S and GOODS-N, respectively, for a total 
area of ~0.9 square degrees. Save for few galaxies, we obtain reliable optical identifications and redshifts, providing us a rich and 
robust dataset for our luminosity function determination. The final combined reliable sample includes 3029 sources, the fraction of 
photometric redshifts being 72% over all redshifts, and almost all at z > 1.5. Based on the multi-wavelength information available in 
these areas, we constrain the LFs at 8, 12, 15 and 24 pm. We also extrapolate total IR luminosities from our best-fit to the observed 
SEDs of each source, and use this to derive the bolometric (8- 1000pm) LF and comoving volume emissivity up to z ~ 2.5. 
Results. In the redshift interval < z < 1, the bolometric IR luminosity density evolves as (1 + z ) 3 8±04 . Although more uncertain 
at higher-z, our results show a flattening of the IR luminosity density at z > 1 • The mean redshift of the peak in the source number 
density shifts with luminosity: the brighest IR galaxies appear to be forming stars earlier in cosmic time (z > 1.5), while the less 
luminous ones keep doing it at more recent epochs (z ~ 1 for L 1R < 10 U L Q X in general agreement with similar data in the literature. 
Conclusions. Our results suggest a rapid increase of the galaxy IR comoving volume emissivity back to 1 ~ 1 and a constant average 
emissivity at z > 1 . We also seem to find a difference in the evolution rate of the source number densities as a function of luminosity, 
a downsizing evolutionary pattern similar to that reported from other samples of cosmic sources. 

Key words. 



1. Introduction 
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A remarkable property of the infrared (IR) selected extragalactic 
sources is their very high rates of evolution with redshift, exceed- 
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ing those measured for galaxies at other wavelengths. Indeed, 
based on data taken by the Infrared Astronomical Satellite 
(IRAS), strong evolution was detected already within the lim- 
ited < z < 0.2 redshift range at far-IR wavelengths (60 pm) 
by Hacking et al. (1987), Franceschini et al. (1988), Saunders et 
al. (1990). The IRAS all-sky coverage also prompted Fang et al. 
(1998) and Shupe et al. (1998) to established the local bench- 
marks at mid-IR wavelengths (12 and 25 pm). The Infrared 
Space Observatory (IS O) allowed us for the first time to perform 
sensitive surveys of distant IR sources in the mid- and far-IR, up 
to z ~ 1 (Elbaz et al. 1999; Puget et al. 1999), showing that dust- 
enshrouded starbursts have undergone strong evolution both in 
luminosity and in density up to that redshift (e.g. Franceschini 
et al. 2001, Chary & Elbaz 2001, Elbaz et al. 2002, Pozzi et al. 
2004). 

The ISO results have been further extended in redshift, to 
Z > 1, by the Spitzer Space Telescope (Papovich et al. 2004; 
Marleau et al. 2004; Dole et al. 2004, Lagache et al. 2004), 
operating between 3.6 and 160 pm with improved sensitivity 
and spatial resolution compared to previous IR observatories. 
Le Floc'h et al. (2005), in particular, studied the evolution of 
IR-bright sources up to z ~ 1 using a sample of mid-infrared 
(24 pm) sources with complete redshift information and derived 
the rest-frame 15 pm and total IR LFs. Strong evolution of the 
IR-selected population with look-back time both in luminosity 
and in density was found in agreement with previous results. 
Similarly, Perez-Gonzalez et al. (2005) have fitted the rest-frame 
12 pm in various redshift bins in the range < z < 3, based on 
photometric redshift determinations, therefore constraining the 
evolution of IR-bright star-forming galaxies. Caputi et al. (2007) 
present the rest-frame 8 pm LF of star-forming galaxies in two 
bins at z ~ 1 and z ~ 2. Babbedge et al. (2006), based on an 
early reduction of the SWIRE EN1 field dataset and using pho- 
tometric redshifts based on relatively shallow optical imaging, 
computed the mid-IR (8 and 24pm) galaxy (AGN1) LF up to 
z ~ 2 (z ~ 4; albeit with uncertainties strongly increasing with 
redshift). Spitzer (Werner et al. 2004) studies have also allowed 
accurate determinations of the local (or low-redshift) LFs at 8 
and 24 pm (Huang et al. 2007; and Marleau et al. 2007; Vaccari 
et al. 2009, in preparation), respectively, essential for compari- 
son with properties of high-z sources. 

All these studies indicated a strong evolution of the IR popu- 
lation up to z ~ 1 (with a comoving luminosity density produced 
by luminous IR galaxies more than 10 times larger at z ~1 than 
in the local Universe). The evolution flattens at higher redshifts, 
up to the highest currently probed by Spitzer, z ~ 3. Moreover, 
luminous and ultra-luminous IR sources, LIRGs and ULIRGs, 
become the dominant population contributing to the comoving 
infrared energy density beyond z ~0.5 and make up 70% of the 
star-forming activity at z~l (Le Floc'h et al. 2005). 

With the aim of achieving an improved knowledge of the 
evolutionary properties of 5p/fzer-selected sources, we exploit 
in this paper the combination of data from the GOODS and 
VVDS-SWIRE multi-wavelength surveys to determine mid-IR 
and bolometric luminosity functions (and thus estimate the SFR 
density) over a wide redshift interval, < z < 2.5. The combina- 
tion of data from three distinct sky areas in our analysis allows us 
to strongly reduce the cosmic variance effects. The coverage of 
the luminosity -redshift plane is also improved with the combina- 
tion of surveys with two limiting fluxes, the deeper GOODS sur- 
veys (80 piy limit) and the shallower wider-area VVDS-SWIRE 
(400 piy limit) survey, with similar number of sources in the 
two flux regimes. The redshift information has been maximized, 
including a large number of new spectroscopic redshifts from 



public databases in the GOODS fields, and an highly optimized 
photometric redshift analysis in the wide-area VVDS-SWIRE 
field, taking advantage of the deep multi-wavelength photome- 
try available and of the VVDS-based photometric redshift tools. 
The source statistics (more than 3000 objects in the combined 
sample) is rich enough for a tight redshift binning in the LF de- 
termination. Finally, the extremely rich suite of complementary 
data from the UV to the far-IR available in the three fields allows 
us to obtain a robust characterization of source spectra in order 
to compute reliable K-corrections and spectral extrapolations to 
be applied to our LF determination. 

At the same time, the large complement of multi-band pho- 
tometric data available for all our sample sources prompts us to 
compute reliable bolometric corrections from the rest-frame 24 
pm to the total (8-1000//m) luminosities. These correction fac- 
tors are particularly well established on consideration that the 
24 pm flux is an excellent proxy to the bolometric flux for active 
galaxies at high-redshifts (for high-z sources longer wavelength 
fluxes are not usually available), see e.g. Caputi et al. (2007), 
Bavouzet et al. (2008). For these reasons we believe that our 
constraints on the evolutionary comoving luminosity density and 
star-formation rate that will be discussed in the paper are among 
the most solid available at the present time. 

The paper is structured as follows. In Section [2] we intro- 
duce the VVDS-SWIRE 24 pm dataset which was first employed 
here. In Section [3] we describe the multi-wavelength identifica- 
tion and redshift determination process for 24 pm sources in the 
GOODS and SWIRE fields. In Section we explain the pro- 
cedure adopted to derive the monochromatic mid-IR rest-frame 
luminosities. Section 5 is devoted to the description of the com- 
putation of the LFs with the 1/V max technique and to the pre- 
sentation of our results as compared to already published data in 
various IR bands. The bolometric LF is also presented. Finally, 
in Section 6 we discuss our results about the evolution of the 
bolometric LF and of the IR luminosity density with redshift. 
Moreover, we compare our results with model predictions. A 
summary of the paper is presented in Section 7. 

We adopt throughout a cosmology with H^-IQ km s" 1 
Mpc~', £2m=0.3 and £2a=0.7. We indicate with the symbol L24 
the luminosity at 24 pm in erg/s (vL v , and similarly for other 
wavelengths). 

2. Spitzer Observations of the VVDS-SWIRE Field 

The 02 hour field portion of the VIMOS VLT Deep Survey 
(VVDS, Le Fevre et al. 2004 and 2005) lies inside the XMM- 
LSS area of the Spitzer Wide-area InfraRed Extragalactic survey 
(SWIRE, Lonsdale et al. 2003, 2006). The total area jointly cov- 
ered on the sky by the two surveys (hereafter, the VVDS-SWIRE 
area) is ~0.85 square degrees. 

The Spitzer observations of the six SWIRE fields were car- 
ried out between December 2003 and December 2004. The 24 
pm observations were executed using the MIPS Scan Map AOT 
with a medium scan rate. Two passes separated by half a field- 
of-view were carried out to allow for removal of cosmic rays and 
transient sources (e.g. asteroids), providing a total exposure time 
of at least 160 s per point, while overlap between rotated scans 
yielded a higher coverage in portions of each map. 

The raw data were reduced using the Spitzer Science Center 
(SSC) standard pipeline. The data processing started from the 
Basic Calibrated Data (BCD) and median filtering was employed 
to even out variations in the local background. While such filter- 
ing optimizes faint source detection, it might also cause a cer- 
tain loss of information about the extended background level in 
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the field. However, this only affects extended or bright (S 24 > 
0.2 Jy) sources, and is thus not an issue for our sample, which 
is mostly made up by point-like faint sources. A slightly more 
sophisticated scheme than in other SWIRE fields had to be put 
in place to remove latent artifacts in the XMM-LSS field caused 
by the proximity of the bright Mira star. Corrected BCDs were 
coadded into large mosaics using the SSC's MOPEX package. 
Source extraction and photometry is described in Section [3~T1 

Shupe et al. (2008) have recently presented the analysis of 
the 24 /mi source counts in the SWIRE fields, with a (differ- 
ential) completeness of 75% at S(2Apm) = 400 /dy. In the 
present work, however, we employed a deeper source catalog 
produced as part of the SWIRE Final Catalog (Surace et al., in 
prep). This catalog is based on the same SWIRE dataset used by 
Shupe et al. (2008) but extracts sources using APEX (a software 
tool provided by SSC) rather than SExtractor (formerly used by 
Shupe et al. 2008) down to a lower SNR. 10.5 arcsec diame- 
ter for point sources, and the Rron flux for extended sources. 
Aperture corrections were computed comparing values obtained 
with smaller apertures and a 30.6 arcsec diameter aperture. Both 
aperture and Rron fluxes were then corrected by the 1.15 factor 
recommended by SSC to take into account the fraction of light 
outside the 30.6 arcsec aperture. We adopted APEX default flux 
measurements based on PRF fitting, as recommended by SSC 
and detailed by the MIPS Data Handbook (SSC, 2007, Version 
3.3.1), and calibrated measured fluxes based on Engelbracht et 
al. (2007). Following the MIPS Data Handbook, we then multi- 
plied APEX fluxes by 1.15 as an aperture correction, to take into 
account the fraction of the PSF that is scattered at large radii, 
and then divided them by 0.961 an a color correction appro- 
priate for an SED with a constant vS v . The differential com- 
pleteness of this deeper catalog is estimated to be better than 
90% at S(24pm) = 400 /Jy in the XMM-LSS area, while its 
reliability at this lower SNR is ensured by the presence of a 
close optical and IRAC counterpart (see Section l3~TT >. We there- 
fore adopted this relatively shallow, but highly reliable, 24 //m 
flux limit in the VVDS-SWIRE region, as a wide-area comple- 
ment to the deeper GOODS samples described in the following 
Section. Our 24 //m reference catalog contains 1494 sources with 
S(24pm) >400//Jy. 

3. Multiwavelength Identification and Analysis of 
the Spitzer MIPS 24 //m Sources 

Our derivation of MIR and bolometric luminosity functions in 
the present paper is based on a Spitzer MIPS 24 /mi sample se- 
lection. We detail in this Section various aspects of the complex 
procedures for source identification in our three reference fields, 
the VVDS/SWIRE, GOODS-S and GOODS-N fields, covering 
a total area of ~0.9 square degrees. 



3.1. VVDS-SWIRE 

We have obtained an observational SED for each one of the 
24 //m sources in the VVDS-SWIRE region by exploiting the 
rich multiwavelength information available from the VVDS 
databas43- The latter includes accurate photometry in the opti- 
cal/NIR spectral range (see description in Arnouts et al. 2007 
and de la Torre et al. 2007), including the deep multicolor pho- 
tometry (u*g'r'i'z') from the Canada-France-Hawaii Telescope 
Legacy Survey (CFHTLS-D1) and the VIMOS VLT Deep 
Survey (VVDS), consisting of deep photometry (Le Fevre et 
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Fig. 1. Fraction of S 24 > 80 //Jy sources with an identified op- 
tical counterpart as a function of the B, V, i and z magnitudes. 
The results are presented both for the separate and combined 
GOODS-S and GOODS-N fields (see the legend in the figure 
for details). The vertical lines mark the 5-cr limits of the optical 
bands. 



al. 2005) in the B,V,R,I,J,Kbands and VIMOS spectroscopy. 
J and K data are also available from the UKIDSS Ultra Deep 
Survey (Lawrence et al. 2007) based on the DR1 release (Warren 
et al. 2007). The SWIRE IRAC photometry is based on the band- 
merged catalog including 3.6, 4.5, 5.8, and 8.0 /an passbands 
(Surace et al. 2005), with a typical 5cr depth of 5.0, 9.0, 43, and 
40 /(Jy respectively. As in Arnouts et al. (2007), we used the flux 
measurements derived in 3 arcsec apertures for faint sources, 
while we adopted adaptive apertures (as for Rron magnitudes 
in Bertin & Arnouts 1996) for bright sources (;7?ab(3.6) < 19.5). 

Photometric redshifts computed by Arnouts et al. (2007) are 
also available through the VVDS database (see also McCracken 
et al. 2003, Radovich et al. 2004, Iovino et al. 2005). These have 
been estimated using the x 1 fitting algorithm Le Phar^l and cal- 
ibrated with the VVDS first epoch spectroscopic redshifts (pro- 
viding ~ 1500 secure spectra for sources with niABO .6fim) < 
21.5) in the same dataset as described in Ilbert et al. (2006), 
including in addition the 3.6 /mi and 4.5 /mi infrared photom- 
etry from SWIRE (Lonsdale et al. 2003). The accuracy achieved 
in the photometric redshift determination is of cr[Az/(l + z)] ~ 
0.031 with no systematic shift. They have been measured be- 
tween < z < 2 (knowing that the high-z tail at z > 1 .4 could 
not be verified with spectroscopic data at that time (Arnouts et 
al. 2007). 

The cross-correlation of the SWIRE 24 /mi catalog with 
the VVDS multiwavelength database was carried out using a 
nearest-neighbor association with a 2 arcsec search radius, since 
at this bright flux levels (S [24pm] > 400 //Jy) confusion is not 
a major issue as it is in the GOODS case. The entire sample 
of 1494 sources detected at 24 /mi over the VVDS-SWIRE re- 
gion was thus assigned an optical/NIR counterpart. For 69 of 
these, however, a photometric redshift (neither spectroscopic or 
photometric) was not available through the VVDS database, and 
we computed it ourselves using Hyperz (Bolzonella et al. 2000) 
and the VVDS-SWIRE photometry as done in Franceschini et 
al. (2006) and in the GOODS case. Given that none of these 
69 sources has a spectroscopic redshift, we are not ableto pro- 
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vide an estimate of the error on our computed photometric red- 
shift. Clearly, these redshifts are much more uncertain because 
their identification was more uncertain. They should be con- 
sidered with some care, but luckily their fraction (4% of the 
VVDS/SWIRE sample, 2% of the total) does not effect the bulk 
of our results. 

All the 1494 sources detected at 24 pm falling within the 
Optic al/NIR/IR AC field have a robust counterpart and redshift. 
The final numbers for the spectroscopic and photometric redshift 
sample are 137 and 1357 (9% and 91%), respectively. 

3.2. GOODS-S Field 

The southern GOODS field, or GOODS-S, is arguably the sky 
area best-studied to date, with deep imaging available from the 
X-rays to the radio and very extensive spectroscopic follow-up. 
Spitzer, in particular, has very deeply surveyed the GOODS-S 
field between 3.6 and 8.0 pm using IRAC and over the range 
24- 1 60 pm using MIPS . 

In this work we made use of the Spitzer MIPS 24 pm 
GOODS-S data products publicly released by the Spitzer 
Science Center as part of the GOODS Spitzer Legacy Data 
Products Third Data Release (DR3) of November 2005 Q. This 
public dataset includes a calibrated map and a catalog of 24 pm 
sources reliable and complete down to S (24pm) - 80 pJy (80% 
complete, as reported for example by Le Floc'h et al. 2005). 
The photometry is based on a PSF fitting algorithm, where the 
SExtractor positions of IRAC 3.6 pm sources are used for input 
to the MIPS source extraction process. The MIPS 24 pm PSF 
was generated from isolated sources in the image, and renormal- 
ized based on the aperture corrections published in the MIPS 
Data Handbook (v2.1, Section 3.7.5, table 3.12). 

In order to build a multiwavelength photometric and spec- 
troscopic catalog for most of the 24 pm sources, we employed a 
number of publicly available catalogs. In the first place, we used 
GOODS-MUSIC (GOODS Multiwavelength Southern Infrared 
Catalog, Grazian et al. 2006, hereafter MUSIC), a 14-band 
multicolor catalog extracted from surveys conducted over the 
GOODS-S region, and summarized in Giavalisco et al. (2004). 
The MUSIC catalog is characterized by excellent reliability and 
completeness, and includes a combination of photometric mea- 
surements that extends from U to 8.0 pm, including £/-band 
data from the ESO 2.2m WFI and VIMOS, the F435W, F606W, 
F775W, and F850LP (BViz) ACS images, the JHK, VLT data, 
and the Spitzer data provided by the IRAC instrument at 3.6, 4.5, 
5.8, and 8.0 pm. 

The MUSIC catalog is jointly selected in the z and K s bands, 
meaning that a z primary selection is performed, followed by a 
K s secondary additional selection. The latter, in any case, adds a 
limited number of sources to the catalog, so that for most prac- 
tical purposes it helps to think of the sample as simply being 
z-selected. Based on this joint z/^-selection and on the posi- 
tions thus determined, fluxes for each source are provided by the 
MUSIC catalog in all available bands, from the UV to the IRAC 
3.6-8.0 pm channels. In the few cases when an IRAC source in 
the field appeared to be undetected in z or K s , we referred to the 
IRAC-selected catalog by Franceschini et al. (2006). Crucially, 
the MUSIC catalog also includes spectroscopic (when available, 
especially from Vanzella et al. 2005 and 2006, and Le Fevre et 
al. 2004) and photometric redshifts for most of the sources. As 
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Fig. 2. Observed redshift distributions in our surveys in 0.2 
redshift bins. Top panel: the separate distributions in the two 
GOODS fields. Middle panel: distributions in the GOODS and 
VVDS-SWIRE fields. Lower panel: continuous line: the redshift 
distribution for the total sample; dot-dashed line: spectroscopic 
redshifts; dashed line: photometric redshifts. 



shown by Grazian et al. (2006), the quality of the resulting pho- 
tometric redshifts is excellent, with an rms scatter in z/(l + z) 
of 0.06 (the same as we achieved in VVDS/SWIRE) and no sys- 
tematic offset over the whole redshift range < z < 6. 

The full GOODS-S 24 pm catalog includes 948 sources over 
an area of about 250 arcmin 2 and is fully reliable and complete 
to S(24pm) = 80 pJy. The MUSIC catalog provides UBVizJK, 
+ IRAC photometric coverage over a 143.2 arcmin 2 area, but, in 
order to maximize the size of our 24 pm sample, we extended 
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Fig. 3. The separate redshift distributions in the two GOODS 
fields are compared with the distribution from Le Floc'h et al. 
(2009), that is based on a large 24 pm sample covering an area 
of 2 square degrees in the COSMOS field at the same flux level 
of our GOODS survey, 5(24) > 80 yuJy. The error bars of Le 
Floc'h et al. account for the effects of cosmic variance. 
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Fig. 4. Combined redshift distribution in GOODS-S & GOODS- 
N fields. This is shown at different flux levels and compared with 
the model of Franceschini et al. (2009). 



this by including areas where either J or K coverage was not 
available (albeit with both UBViz and IRAQ. The final effective 
area of our 24 //m catalog is thus of 168 arcmin 2 , and the fi- 
nal sample totals 666 sources. For each of these we searched for 
counterparts in the MUSIC catalog. Given the well-established 
one-to-one relation between 3.6 and 24 fim sources, and the 
power demonstrated by the deep IRAC 3.6 yum observations to 
identify even the deepest 24 fim sources (see e.g. Rodighiero et 
al. 2006), we have based our 24 fim identification process on the 
MUSIC 3.6 fim catalogue. 

The likelihood and the reliability of the associations were 
computed using the likelihood ratio (hereafter simply likelihood) 
technique described by Sutherland & Saunders (1992), whose 
practical implementation is detailed in Vaccari et al. (2009, in 
preparation). According to this, the likelihood is defined as the 
ratio between the probability that an IRAC 3.6 pm source is the 
real counterpart of a MIPS 24 yum source and the probability 
that the association is actually with a nearby background source, 



under the assumption that magnitudes and positions are not cor- 
related (i.e. no constraint is set on the source/association flux 
ratio). 

For each 24 yum source, we considered all 3.6 /jm sources 
within a search radius of 5 arcsec as potential counterparts. This 
is approximately 3 times the quadratically summed rms radial 
positional errors of MIPS and IRAC sources, estimated to be 1 .5 
and 0.75 arcsec rms respectively. These values are a bit larger 
than generally assumed but are believed to be reasonable given 
that at the flux limits we are considering the formal positional er- 
rors do not account for confusion effects. The likelihood and re- 
liability values are thus determined by the source number counts 
of the association catalogue and the source-target relative posi- 
tions on the sky. 

Likelihood and reliability were then computed for all poten- 
tial counterparts. Out of 666 sources, 499 (75%) were assigned 
a counterpart and a redshift based on the MUSIC catalog which 
was both highly likely and reliable (>99%), whereas the remain- 
ing 167 (25%) were flagged for visual checks (i.e. having lower 
likelihood and/or reliability, or no potential counterpart within 
the 5 arcsec search radius). We also tried additional likelihood 
and reliability analyses with reference to the z- and K s entries 
of the MUSIC catalog, and verified that the latter association did 
not produce any substantial changes and therefore did not pursue 
this any further. 

The visual checks helped us to solve or at least explain most 
of the ambiguities. Occasionally, the MUSIC catalog turned out 
to be badly hampered by confusion, or its z/K s joint selection 
missed some clearly detected IRAC 3.6 /um source. In most 
cases, however, visual checks provided a straightforward iden- 
tification with a given source, or with a close pair of sources at 
very similar redshifts and thus likely in the process of merging. 
The remaining cases, totalling less than 10% of the sample and 
generally consisting of the most confused, disturbed or faintest 
sources, proved somewhat of a challenge and required querying 
a number of different catalogs, as well as implementing our own 
photometric redshift solution. More specifically, we combined in 
these cases information available through the COMB017 cata- 
log (Wolf et al. 2004) and our own IRAC 3.6 yum selected catalog 
(Franceschini et al. 2006) to derive the full optic al/NIR/IRAC 
SED of these 24 yum sources, and determine a reliable redshift 
running an updated version of Hyperz, following the approach 
detailed by Franceschini et al. (2006). 

Out of the 167 (25%) ambiguous cases above, 59 (9%) were 
assigned a robust counterpart and redshift from the MUSIC cat- 
alog following visual checks, 56 (8%) from the COMB017 cat- 
alog, and 32 (5%) from the Franceschini et al. (2006) catalog, 
leaving only 20 sources, or just 3% of the original sample, with- 
out a redshift. All of these extra redshifts were verified and oc- 
casionally corrected when reported as unreliable in the original 
catalogs by running Hyperz on the full SEDs. 

During the progress of this work, when the identification 
and SED analysis of GOODS-S sources had already been com- 
pleted, a number of additional spectroscopic redshifts were made 
available (Ravikumar et al. 2007, Vanzella et al. 2008, Popesso 
et al. 2009), and Wuyts et al. 2008 compiled them into their 
FIREWORKS catalog. We have then adopted the FIREWORKS 
spectroscopic redshifts for the optical/NIR/IRAC sources which 
we had previously identified as robust counterparts to 24 yum 
sources. We have used a very stringent 1 arcsec search radius 
for the association of our MIPS/IRAC/MUSIC sources with 
the FIREWORKS catalogue, hence obtaining 165 extra spectro- 
scopic redshifts for sources for which only photometric redshifts 
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were previously available, greatly adding to the overall quality 
of our final sample. 

Eventually, out of the 666 sources detected at 24 pm within 
the ACS+IRAC field, 646 (97%) sources were assigned a robust 
counterpart andredshift, with only 20 (3%) unidentified sources. 
The final numbers for spectroscopic and photometric redshifts 
are 415 and 231 (64% and 36%), respectively. 

3.3. GOODS-N Field 

The northern GOODS field, or GOODS-N, was observed by 
HST ACS and by Spitzer IRAC and MIPS at virtually the same 
depths and over a similar area as GOODS-S, thus providing a 
useful comparison to evaluate cosmic variance and to improve 
the statistics at faint flux levels. In the literature, however, the 
field was so far given less attention than its southern counter- 
part, and while data have been publicly available for some time, 
no multi-wavelength public catalogs have been released to date. 
In this work we have made use for our primary selection of 
the Spitzer MIPS 24 pm GOODS-N data products publicly re- 
leased by the Spitzer Science Center as part of the GOODS 
Spitzer Legacy Data Products Interim Data Release 1 (IDR1) 
of February 2004. In order to carry out a multi-wavelength iden- 
tification and SED fitting of the 24 pm sources as done in the 
previous subsection, we used an original catalog containing full 
photometric information (U from MOSAIC at KPNO-4-m, BViz 
from ACS, JHK S from FLAMINGOS at KPNO 4m, IRAC chan- 
nels), as well as a detailed compilation of spectroscopic and pho- 
tometric redshifts (Mancini et al. 2009). After removing 5<x out- 
liers of the best-fit solution, Mancini et al. found a combined 
mean offset of (z spec -z p hot)/(l +z sp ec) = 0.004 and arms scatter 

Of CT [(z spec - Zp hot )/(l + Z S pec)] = 0.09. 

Such a comparison catalog was selected in the IRAC 4.5 pm 
band, which provides image quality and depth similar to IRAC 
3.6 pm. This pure IRAC selection reduces the source confusion 
problems with respect to a z- or Zfj-selected comparison sam- 



ple: on one side, most of the disturbed z/K s source pairs appear 
as unresolved in IRAC images, on the other because some ex- 
tremely faint z/K s sources are more robustly detected at IRAC 
wavelengths, both factors contributing to an easier identification 
process. 

We thus carried out a likelihood and reliability analysis of 
associations between the 24 pm and 4.5 pm sources in the same 
way as it was done in the previous Section, but the number 
of ambiguous identifications was now extremely limited. In the 
process, we also verified that using the IRAC 3.6 pm and IRAC 
4.5 pm comparison samples for identification purposes produced 
virtually identical results. 

Of the 904 24pm sources of the primary selection falling 
within the ACS+IRAC field, only 26 (3%) were flagged for vi- 
sual checks by the likelihood and reliability criteria outlined in 
the previous Section, and 889 (or 98.3%) sources found a robust 
counterpart and redshift after the visual checks, with only 15 
(1.7%) unidentified sources. The final numbers of spectroscopic 
and photometric redshifts are 290 and 599 (32.5% and 67.5%), 
respectively. 

In Figure Q] we report the fraction of S24 > 80 p]y sources 
with an identified optical counterpart as a function of the B, V, i 
and z magnitudes. The results are presented both for the separate 
and combined GOODS fields. 



4. Redshift distribution and source number counts 

To better characterize the statistical properties of our samples, 
we report in Figure [2] the redshift distributions for the 24 pm 
sources. We compare there the separate distributions in GOODS- 
N and GOODS-S and the GOODS with WDS-SWIRE source 
samples at the respective flux limits of S(24pm) > 80 pJy 
and S{24pm) > 400 pJy. As apparent in the figure, our final 
source sample not only guarantees excellent statistics over the 
< z < 1.5 redshift interval, but also includes a substantial sta- 
tistical coverage of the high redshift population, up to z ~ 3, 
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which is an absolute novelty of the Spitzer mission for cosmol- 
ogy compared to the previous ISO surveys. In Figure [3] we also 
report the separate distributions in the two GOODS fields com- 
pared with the distribution from Le Floc'h et al. (2009) that is 
based on a large 24 pm sample covering an area of 2 square de- 
grees in the COSMOS field at the same flux level of our GOODS 
survey (i.e. 80 pJy). The three distributions have been normal- 
ized to the same sky area, 1 square degree. The COSMOS dis- 
tribution has been renormalized to the GOODS area. The error 
bars of Le Floc'h et al. account for the effects of cosmic variance, 
showing that the combined use of the two GOODS sky areas is 
very useful to knock down any local inhomogeneity. In Figure|4] 
we show again the combined redshift distribution in GOODS-S 
& GOODS-N fields, this time splitted at different flux levels and 
compared with the model of Franceschini et al. (2009). 

The impact of cosmic variance in our work can be also veri- 
fied by looking at the extragalactic source number counts in the 
various sub-samples we are using. This is presented in Figure 
[5] where we show the 24 pm differential number counts in the 
GOODS-S, GOODS-N & VVDS-SWIRE fields, separately, and 
compared with predictions from the Franceschini et al. (2009) 
model, and with previous measurements by Chary et al (2004), 
Papovich et al. (2004) and Shupe et al. (2008). We observe a 
general agreement among the various surveys, with the excep- 
tion of the GOODS-S that below S (24pm) < 0.5 mJy lies below 
the other data-set. A similar trend can be also observed in the 



redshift distribution (Figure 2, top panel) and later in the 24 pm 
luminosity function (Figure 7) and should be really attributed to 
the cosmic variance effects. 



5. Infrared luminosities from multiwavelength SED 
fitting 

The rich multiwavelength optical-to-IR dataset available in the 
GOODS and VVDS-SWIRE fields is ideal to perform a spec- 
tral fitting procedure on the whole observational SED of each 
source in our sample. This is needed for the calculation of the 
K-corrections to estimate the rest-frame luminosities at different 
frequencies and the IR bolometric corrections. Since we already 
know the redshift of each object from the spectroscopic or pho- 
tometric measurement, we used Hyperzspec (Bolzonella, private 
communication), an adapted version of Hyperz performing SED 
fitting at a fixed redshift, over the whole broad-band photomet- 
ric set available to us. In this fitting procedure, we compared the 
optical-to-IR SEDs of 24 pm sources with a library of template 
SEDs of local objects. The library contains 20 spectra, including 
one elliptical, seven spirals, three starbursts, six AGNs, and three 
composite (starburst + AGN) templates covering the wavelength 
range between 1000 A and 1000 pm (for a detailed description, 
see Polletta et al. 2007 and Franceschini et al. 2005). Hyperzspec 
finds the best-fit SED by minimizing the x 2 derived from com- 
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Fig. 7. L24 rest-frame 24 fim luminosities as a function of redshift for GOODS (red asterisks) and VVDS-SWIRE (blue crosses) 
sources. The thick green solid and the blue dashed lines indicate as a function of redshift the 24 pm rest-frame luminosity corre- 
sponding to an observed 24 pm flux of 0.08 mJy with the template of an SB galaxy and of M82, respectively. The green symbols 
indicate the distribution of Type 1 AGNs as selected from the SED fitting procedure. 



parison of the observed and theoretical SEDs at the redshift of 
each source. 

Given that the limited library we are using can not be ex- 
pected to satisfactorily reproduce the observed SED from the 
UV to the far-IR data of each source, we have forced Hyperz to 
favour the fitting of the IR data (e.g. IRAC and MIPS photomet- 
ric points). We did this by artificially changing the photometric 
errors and providing much smaller values in the Spitzer regime. 
In detail, we adopted the following receipt: for optical data (i.e. 
A < Ipm) we assumed errors of 1 mag, for near-IR data ( J,H 
and K bands) errors of 0. 15 mag and for Spitzer IRAC and MIPS 
data errors of 0. 1 mag. This combination of weights provided the 
best set of fit to the sources, generally always providing a good 
description of the mid- to far-IR observed shape of each SED. 
For the purposes of this work, a poor fit of the optical side of the 
spectrum for some sources would not be critical, as we are inter- 
ested in deriving luminosities at rest-frame wavelengths above 8 
pm. However, the inclusion of the whole SED in the fitting pro- 
cedure allow us to fully exploit the photometric information and 
better characterize the K-correction. 

This procedure provided us with a best-fit SED for each 
galaxy (the few stars where removed a priori on the basis of their 
optical morphology and/or IR colors) and with a rough classifi- 
cation between star-formation and AGN-dominated activities as 
the main source of the far-IR emission. As an example, we re- 
port in Figure [6] the best-fits that we obtained for four sources 
representative of the main spectral classes considered. One of 
the best-fit is with a typical spiral (Sb) spectrum, one with a star- 



burst galaxy (M82), a source whose IR emission is contributed 
both by star-formation and nuclear activity (Markarian 231), and 
finally a type-1 quasar (QSO). 



In particular, we identified as type-1 AGN those sources clas- 
sified by the automatic tool as type-1 quasars or Seyfert-1, cor- 
responding to five SED templates in the Polletta et al. (2007) 
set. These are characterized by monotonically rising spectra (e.g. 
Alonso-Herrero et al. 2005), similar to the QSO spectrum in Fig. 
[6] We found that this class of objects includes 320 sources, or 
~10% of our complete sample. However, this fraction only con- 
cerns the bright end of our IR LFs, and does not affect our results 
about their evolution. For this reason, we excluded them from 
any statistical analysis in the following part of the paper, like in 
the computation of the LFs. 



The identification of type-2 AGNs or other active sources 
is much more uncertain, given the complexity of their spectra 
in the presence of mixed contributions from star-formation and 
nuclear activity. The lack of photometric information in the 10- 
20 pm region of the SED prevented us to be more specific in 
this classification process in the (majority) case of a combined 
AGN/starburst dust-obscured activity (see also Gruppioni et al. 
2008). We defer to a future work for a more detailed analysis 
of the AGN luminosity functions based on the present sample, 
using color-color diagnostics and X-ray information. 
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Fig. 8. Rest-frame 24 pm luminosity function for all galaxies in the GOODS-North and GOODS-South separated, as computed with 
the 1 IV„ mx . As a reference, we report the local luminosity function at 25 pm (dotted green line) by Shupe et al. (1998). The red open 
circles represent the original estimates of the GOODS-N 24mu LF values, together with their poissonian error bars. The red upper 
and lower solid lines represent the range of values derived with 100 iterations by allowing a change in photo-z and K-correction. 



For each source in our combined sample, we calculated the 
rest-frame NIR and MIR luminosity densities as: 

And 2 , 

L v (v rest ) = S v (v obs ) (1) 

1 +z 

where di is the luminosity distance for a given redshift in our 
adopted cosmology, S v the flux density in erg crrT 2 s _1 Hz -1 , 
and v bs and v rest are the observed and rest-frame frequencies, 
respectively, where v rest — (1 + z)v b s - Throughout the paper, we 
present IR luminosities in terms of the bolometric luminosity of 
the Sun, L = 3.83 x 10 33 erg s _1 . 

The distribution of the 24 pm luminosities as a function of 
redshift for all the sources in the GOODS (red symbols) and 
VVDS-SWIRE (blue symbols) samples are reported in Figure[7] 
(type 1 AGNs are marked as green symbols). This figure shows 
that the combination of the two samples offers a fair coverage 
in luminosity up to z ~ 2.5. To show the completeness of our 
survey we also report as thick green solid and blue dashed lines, 
as a function of redshift, the 24 pm rest-frame luminosity corre- 
sponding to an observed 24 pm flux of 0.08 mJy with the tem- 
plate of an SB galaxy and of M82, respectively. As can be now 
clearly seen in Figure 4, the bimodal distribution in the 24 /mi lu- 
minosity as a function of redshiftis an effect of the K-correction 
(quiescent spiral versus starburst galaxy). 

6. Luminosity Functions 

A variety of rest-frame wavelengths has been adopted by differ- 
ent authors to compute from Spitzer data the redshift evolution of 



the MIR luminosity functions and to compare them to the local 
ones. Caputi et al. (2007) have recently presented the rest-frame 
8 pm LF at z ~ 1 and z ~ 2, Perez-Gonzalez et al. (2005) derived 
the rest-frame 12 pm LF up to z ~ 2.5, Le Floc'h et al. (2006) 
obtained the 15 pm LF up to z ~ 1.2 and Babbedge et al. (2006) 
finally computed the LF at both 8 and 24 pm up to z ~2, albeit 
with uncertainties substantially increasing with redshifts arising 
from their relying on shallow Spitzer and optical photometry and 
on purely photometric redshifts. Some authors also tried to per- 
form the computation of the evolutionary IR [8-1000^m] LF. 

Our largely improved sample allows us homogeneous selec- 
tion and treatment of a large database, including data from the 
wide VVDS-SWIRE field, which are essential for a proper sam- 
pling of the high-luminosity end of the LF and for keeping under 
control the cosmic variance problems of the deeper samples. In 
turn, the GOODS samples (smaller but a factor 5 deeper in flux, 
S (24) > 80 fi}y) have already shown their power to constrain 
the evolution of the LF at lower luminosity levels and higher 
redshifts (e.g. Caputi et al. 2007, Le Floc'h et al 2005). 

6.1. LF computation 

The LFs at the rest-frame wavelengths of 8, 12, 15 and 24 pm 
from our combined GOODS+VVDS-SWIRE sample have been 
computed using the 1 /V max method (Schmidt 1968). The advan- 
tage of this technique is that it allows direct computation of the 
LF from the data, with no parametric dependence or model as- 
sumption. We split the sample into redshift bins selected so as 
to offer adequate numbers of galaxies in each. Our adopted red- 
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Fig. 9. Rest-frame LF at 8 /mi. The blue thick lines mark our combined GOODS+VVDS-SWIRE luminosity functions in the various 
redshift bins. In each panel we compare our results with literature data (when available). The green-dotted line represents the local 
LF as computed by Huang et al. (2006). The filled black symbols are from Caputi et al. (2007), while the open pink squares are 
from the SWIRE survey (Babbadge et al. 2006). The vertical dot-dashed lines reported in each redshift bin represent the luminosity 
above which we do not expect any incompleteness. We also report as black-dashed lines the model predictions by Franceschini et 
al. (2009). 



shift binning was different in different wavelengths, in order to 
provide an easier comparison with published data. 

In each redshift bin, the comoving volume available to each 
source is defined as V max = V Zmar - V ZmM . The Zmax value corre- 
sponds to the maximum redshift at which a source would still be 
included in the sample, given the limiting 24 //m flux, while z m ,„ 
is the lower limit of the considered redshift bin. In any case, z max 
should be smaller than the maximum redshift of the bin consid- 
ered. 

As previously mentioned in Section 3, all our 24 /mi cata- 
logs are basically complete down to our adopted limiting fluxes 
(400 and 80 //Jy), and thus we did not apply any completeness 
correction. 

We have checked the consistency, within the statistical errors 
(the error bars include only the Poisson noise), and the effects of 
cosmic variance in our LF by comparing the results obtained in- 
dependently in the separate fields. In particular, the rest-frame 24 
/mi LF calculated from the GOODS-North and in the GOODS- 
South samples are reported in Figure [8] The two LFs look are 
consistent with each other within the error bars (similarly to what 
found by Caputi et al. 2007 for their rest-frame 8 //m LF). 

To smooth the effects of cosmic variance and to consolidate 
the statistics, we provide in the following combined results from 
the three samples. 



In order to better quantify the effects of the photometric red- 
shifts and of the K-correction on our results concerning the esti- 
mate of the luminosity functions, we performed a set of Monte 
Carlo simulations. We used as a test case the GOODS-N sam- 
ple, and checked the effect on the 24 /mi luminosity function, 
that is the one that requires the major extrapolation from the 
data, and so should be the more critical case in our work (sim- 
ilarly to the bolometric LF). We have iterated the computation 
of the 24 fim LF by perturbing each time the photometric red- 
shift of each source (by an amount corresponding to the RMS in 
the photo-z/spec-z relation reported for the GOODS-N, i.e. 0.09) 
Every time we then recomputed the monocromatic luminosity by 
performing, as in the main case, the SED fitting with Hyperz and 
the computation of the V„ mx . Hyperz is left free to choose a new 
best-fit template, so that we take simultaneously into accounting 
for a variation in redshift and K-correction. We included also 
Arp200 in the library of templates to see if the effect of a more 
extreme ULIRG does affect our results. The results of this Monte 
Carlo simulation are also reported in Figure|8]. The red open cir- 
cles represent the original estimates of the GOODS-N 24mu LF 
values, together with their poissonian error bars. The red upper 
and lower solid lines represent the range of values derived with 
100 iterations by allowing a change in photo-z and K-correction. 
The comparison shows that indeed the combined effects of pho- 
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Fig. 10. Rest-frame LF at 12 /jm. The blue thick lines mark our combined GOODS+VVDS-SWIRE luminosity functions in the 
various redshift bins. In each panel we compare our results with literature data (if available). The green-dashed line represents the 
local LF (Rush et al. 1993). The the red-crosses are from Perez-Gonzalez et al. (2005). The vertical dot-dashed lines reported in each 
redshift bin represent the luminosity above which we do not expect any incompleteness. We also report as black-dashed lines the 
model predictions by Franceschini et al. (2009). 



tometric redshift and K-correction on the error bars are larger 
than the simple poissonian values. But looking at the GOODS-S 
data (blue circles) we seem to find that cosmic variance has a 
greater impact than the photo-z uncertainties. Given that we can 
not extensively estimate the cosmic variance in the full survey 
(and that we are using different fields to minimize this problem) 
in the paper we prefer to still use the poissonian error bars to 
derive our conclusions. 



6.2. Results for various IR bands 

In Figures [9] to [12] we present the results of the computation of 
our rest-frame 8, 12, 15 and 24 fim luminosity functions, re- 
spectively. As mentioned above, the choice of the redshift bins 
is different for different wavelengths, depending on the avail- 
able published data which we use for comparison. The redshift 
and luminosity binning is also requested to provide a statisti- 
cally sufficient number of sources. The error bars include only 
the Poisson noise. The vertical dot-dashed lines reported in each 
redshift bin represent the luminosity above which we do not ex- 
pect any incompleteness. This is clearly evident from the data 
themselves, that start to fall down below this level, mainly be- 
cause at fainter luminosities not all galaxy types are observable 
(depending on their SED, see Ilbert et al. 2004 for a detailed 
discussion of this bias). The center of each luminosity bin has 



been computed as the median of the observed luminosities dis- 
tribution in the corresponding bin. The LFs from the total com- 
bined GOODS+VVDS-SWIRE, including all galaxy types, are 
reported in each figure panels with a thick blue line. The val- 
ues of the monochromatic LFs at 8, 12, 15 and 24 fim are also 
reported in Tables[T]|2][3]and|4] respectively. 

Our results in Figs. [9] to [12] are compared with available data 
in the literature, whenever available. At 8 //m (Figure [9]l the 
most significant comparison is with the results of Caputi et al. 
(2007), who computed the rest-frame LF at z ~ 1 and z ~ 2 
from the combined dataset in GOODS-S and GOODS-N. We 
find a fairly good agreement within the errorbars, still with some 
differences, in the two redshift bins, reaching similar luminos- 
ity levels of completeness. However, the addition of the VVDS- 
SWIRE data allow us to extend the LF to higher luminosities. 
This result might confirm the finding of Caputi et al., that dif- 
ferent approaches in the estimation of the /f-correction provide 
consistent results when deriving the monochromatic IR luminos- 
ity functions at high-z. Up to z ~ 1, we also find good consis- 
tency with the results by Babbedge et al. (2006) based on the 
SWIRE survey. However, at higher z there results significantly 
diverge from ours, the LF being much lower. 

At 12 /zm (Figure ITOb a valuable comparison is possible up 
to z ~ 2.5 with the LF derived by Perez-Gonzalez et al. (2005), 
based on a large sample of more than 8000 sources brighter than 
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Fig. 11. Rest-frame LF at 15 fjm. The blue thick lines mark our combined GOODS+VVDS-SWIRE luminosity functions in the 
various redshift bins. In each panel we compare our results with literature data (if available). Open circles represent the local LF 
(Pozzi et al. 2004), red-stars are from Le Floc'h et al. (2005) and the cyan dot-dashed lines are from Magnelli et al. (2009). The 
vertical dot-dashed lines reported in each redshift bin represent the luminosity above which we do not expect any incompleteness. 
We also report as black-dashed lines the model predictions by Franceschini et al. (2009). 



S (24/j.m) = 80 fi]y. Although this work is based on purely photo- 
metric redshifts, we observe a general good agreement with our 
estimates up to z ~ 1.8, for what concerns the bright-end of the 
LFs, while we note that our LFs seem to be lower in the fainter 
luminosity bins computed from our data. In the two higher red- 
shift bins (1.8 < z < 2.2 and 2.2 < z < 2.6) the LF values 
observed in our case are lower, a difference that is likely to at- 
tribute to the combination of the uncertainties in the estimates of 
photometric redshifts and to those in the /f-corrections, and to 
the cosmic variance. 

Similar considerations can be made about the 15 fim results 
(Figure fTTTi. where the major reference comes from Le Floc'h et 
al. (2005) LF estimation up to z ~ 1 . We also report the recent re- 
sults by Magnelli et al. (2009) up to z ~ 1 (we report their best-fit 
double power-law). In this case, our combined GOODS+VVDS- 
SWIRE sample provides consistent results in the common red- 
shift range. A comparison is also shown with the local LF by 
Pozzi et al. (2004). 

Finally, we have attempted to compute for the first time the 
rest-frame 24 fim LF up to z ~ 2.5. The results are presented in 
Figure [12] A comparison can be performed at low redshifts with 
the IRAS 25 yum LF (Shupe et al. 1998), and with the Spitzer re- 
sults from the First Look Survey (FLS) 24 fim team (Marleau et 
al. 2007). As observed at 8 pm, the LF computed by Babbedge 
et al. (2006) is consistent with our results at 24 //m up to z ~ 1 . 



Above, their values seem to underestimate the present compu- 
tation. We also report a recent estimate by Vaccari et al. (2009) 
based on the SWIRE-SDSS database. We note the very good 
agreement of our estimate of the z < 0.3 LF with the local ones. 

6.3. Total infrared bolometric luminosities and LF 

Many efforts have been made by various authors to convert the 
observed MIR luminosities into total infrared (TIR) luminosi- 
ties (e.g. Chary & Elbaz 2001, Elbaz et al. 2002, Takeuchi et al. 
2005, Caputi et al. 2007, see Bavouzet et al. 2008 and references 
therein for a recent review). There is an obvious strong interest 
in attempting to estimate the bolometric luminosities of cosmic 
sources, and the evolution of the bolometric luminosity func- 
tions, because these quantities can be directly scaled to the rates 
of gas processing into stars (Kennicutt 1983; Rowan-Robinson 
et al. 1997), for the stellar-dominated, and of gas nuclear accre- 
tion, for the AGN-dominated objects. Based on studies of the 
spectral energy distributions of local sources, such corrections 
to (our reference) 24 yum fluxes are relatively well established 
(e.g. Caputi et al. 2007; Bavouzet et al. 2008). 

Our adopted procedure in the present paper has been to inte- 
grate for each object the best-fit SED, that we obtained in Sect. 
[5] in the rest-frame [8-1000yum] interval, as a measure of the 
corresponding TIR luminosity. The result of this computation is 
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Fig. 12. Rest-frame LF at 24 /mi. The blue thick lines mark our combined GOODS+VVDS-SWIRE luminosity functions in the 
various redshift bins. In each panel we compare our results with literature data (if available). The dotted-green line represents the 
local LF (Shupe et al. 1998). The black asterisks are the local LF as computed by Marleau et al. (2007), the open pink squares are 
from the SWIRE survey (Babbedge et al. 2006) and the red filled circle shows a recent estimate by Vaccari et al. (2009) based on 
the SWIRE-SDSS database. We also report as black-dashed lines the model predictions by Franceschini et al. (2009). The vertical 
dot-dashed lines reported in each redshift bin represent the luminosity above which we do not expect any incompleteness. 



shown in Figure [13] where we report with different symbols the 
distribution of the TIR luminosities for our GOODS and VVDS- 
SWIRE sources, as a function of redshift. Our approach is sim- 
ilar to others recently reported in the literature (Le Floc'h et al. 
2005, Perez-Gonzalez et al. 2005, Caputi et al. 2007), all us- 
ing similar sets of templates of local IR galaxies to fit the data. 
As in Figure 3, the thick green solid and blue dashed lines in 
Figure 9 indicate, as a function of redshift, the total IR luminos- 
ity corresponding to an observed 24 /mi flux of 0.08 mJy with 
the template of an SB galaxy and of M82, respectively. 

In order to check the consistency of our approach, we com- 
pared our results with those of Takeuchi et al. (2005), Sajina et 
al. (2005), Caputi et al. (2007), and Bavouzet et al. (2008), in 
particular. In their work, Caputi et al. provide an empirical cal- 
ibration of the conversion from vL v (8yum) to the bolometric IR 
luminosity (jJ^X based on a sample of Sp/fzer-selected galax- 
ies, while Takeuchi et al. (2005) and Sajina et al. (2005) obtain 
a relation between vL v (24/«m) and L( s ,. Bavouzet et al. (2008) 
elaborate further the latter by including into the analysis data 
at 8, 24, 70 and 160 /mi, mostly from Spitzer, for their sample 
sources. These relations should then be reliable up to z ~ 2 and 
Z ~ 1 for ULIRGs and LIRGs, respectively. 

In Figure [14] we report a comparison of our observed 
monochromatic versus bolometric relation for all our sample 
sources with the average adopted relations by Caputi et al. (8 
/mi; left panel) and Bavouzet et al. (24 //m; right panel). The var- 



ious datapoint alignments in the plots correspond to the different 
spectral templates that we used to fit our photometric data. 

The solid line in the left panel shows the empirical rela- 
tion provided by Caputi et al. (2007) for the rest-frame 8 /mi 
wavelength. We observe an excellent agreement between our 
average relationship and that by Caputi et al. (2007), over the 
whole luminosity range. The relations provided by Bavouzet et 
al. (2008) are calibrated using longer wavelength MIPS data (70 
and 160 /mi) and should be considered as particularly reliable. 
We report those for the 8 /mi and 24 /zm (dashed and dot-dashed 
lines in both panels of Figure [14l calibrated on local and high- 
redshift galaxies, respectively). At 8 /mi the relation of Bavouzet 
et al. appears to be somewhat flatter than ours, but still consis- 
tent within the data dispersion. However, we have to note that 
when extrapolating the relation at high-z, Bavouzet et al. use the 
stacking technique, and this may get wrong faint IR luminosities, 
indeed their faint total IR luminosities could be overestimated 
compared to our data. 

The comparison of the monochromatic 24 /mi and bolomet- 
ric luminosities in the right panel shows good agreement with 
the relations by Takeuchi et al. (2005) and a slight offset with 
Sajina et al. (2005). On the other hand, the slope of the Bavouzet 
et al. relation appears to be flatter than ours. This difference sug- 
gests a potential underestimate of our TIR luminosities below 
log(L/s)<10.5 L Q in the lowest redshift bin (z < 0.3). However, 
given the flatness of the LF at such faint luminosities, even a sig- 
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Fig. 13. Total IR luminosities as a function of redshift for GOODS (red asterisks) and VVDS-SWIRE (blue crosses) sources. The 
thick green solid and the blue dashed lines indicate as a function of redshift the total IR luminosity corresponding to an observed 24 
/jm flux of 0.08 mJy with the template of an SB galaxy and of M82, respectively. 




Fig. 14. Distribution of total IR luminosities as a function of the corresponding monochromatic rest-frame 8 (left panel) and 24 
jum (right panel) luminosities for GOODS (red asterisks) and VVDS-SWIRE (blue crosses) sources. The solid line in the left 
panel shows the empirical relation provided by Caputi et al. (2007) for the rest-frame 8 /mi wavelength. The relations provided by 
Bavouzet et al. (2008) are reported in both panels as dashed and dot-dashed lines (calibrated on local and high-redshift galaxies, 
respectively). The cyan and green lines show the relations by Takeuchi et al. (2005) and Sajina et al. (2005), respectively. 



nificant underestimate of the true L/r at these luminosities would 
not strongly affect our results. 

By applying the same method as described in Section 16.11 
to compute the monochromatic luminosity functions, we have 
derived the bolometric [8-1000 /mi] LF up to z ~ 2.5. The results 
are presented in Figure Q3] and reported in Table [5] 



As we did for the monochromatic LF discussed in Sect. 16.21 
we have compared our results with those available in the recent 
literature, adapting our choice of redshift bins. Again, we find a 
general consistency with other results at various redshifts, in par- 
ticular with the LF estimates of Le Floc'h et al. (2005), Huynh 
et al. (2007) and Magnelli et al. (2009) up to z ~ 1. The lat- 
ter is based on the Spitzer MIPS 70 /mi data in the GOODS- 
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Fig. 15. Bolometric IR (8-1000/mi) LF. The points connected by blue thick lines mark our combined GOODS+VVDS-SWIRE 
luminosity functions in the various redshift bins (here error bars are plotted in red for clarity). The black solid lines show our 
best-fit using a parameterized function. In each panel we compare our results with literature data (if available). The dashed-green 
line represents the local LF as computed by Sanders et al. (2003). The open-pink diamonds are from Huynh et al. (2007), while 
the open-red stars are from Le Floc'h et al. (2005). The cyan dot-dashed lines are from Magnelli et al. (2009). We also report 
with black-filled circles the results of Caputi et al. (2007) and with black asterisks the data by Chapman et al. (2005). The vertical 
dot-dashed lines reported in each redshift bin represent the luminosity above which we do not expect any incompleteness. 



N field and should be considered as particularly reliable due 
to the moderate bolometric correction. At higher redshifts, we 
still observe a general agreement with the results published by 
Caputi et al. (2007) at z ~ 2 and some consistency also with 
the 2 < z < 2.5 bolometric LF at the extremely high luminos- 
ity values derived by Chapman et al. (2005) from radio-detected 
submillimeter galaxies. 

For visual intuition of the observed LF evolution, we report 
in each panel the local bolometric LF computed by Sanders et 
al. (2003) from the IRAS revised galaxy sample. 

7. Discussion 

7.1. The evolutionary bolometric luminosity functions 

For an easier description of the evolution of the bolometric LF, 
we adopted a routinely used parameterization law (Saunders et 
al. 1990, Pozzi et al. 2004, Le Floc'h et al. 2005, Caputi et al. 
2007) to fit our data: 



<D(L)dlog 10 (L) = 



**(£) " eXp -2^ l0g -( 1 + ^) 



dlog 10 (L) 



(2) 



where, in this case, L is the TIR [8-1000jt/m] luminosity. The 
parameter a correspond to the slope at the faint end. L* is the 
characteristic Zi*. luminosity and O* is the normalization factor. 
The best-fitting parameters for the local LF are: a - 1.2, cr ^ 
0.7, L* 1.77 1O 9 L and <t>* 0.0089 Mpc^. 

The main uncertainties potentially affecting our results are 
those concerning the estimate of the bolometric luminosity from 
multiwavelength data limited to A < 24 /mi, and the lack of 
statistics at the faint end of the LF preventing us from signifi- 
cantly constraining its slope. We have kept the slope of the faint- 
end of the LF fixed to the local observed value (or = 1.2; see 
also Zheng et al. 2006, Caputi et al. 2007) and we have fitted 
the observations in each redshift bin by varying only L* and O*, 
using a x 1 minimization procedure. The slope at the bright end 
has been only slightly changed manually in a couple of redshift 
bins to provide a better fit in the highest luminosity bin. This 
procedure corresponds to a combination of luminosity and den- 
sity evolution. The results of the fitting procedure are presented 
in Figure Q3] (solid black lines). 

From our fitting, we infer that the comoving number density 
of sources, as parameterized by <£*, evolves like (1 + z) in 
the redshift range < z < 1, while the characteristic infrared 
luminosity (L*) evolves as (1 +z) 2 ' 7 - Above z > 1, the degeneracy 
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Fig. 16. Evolution of the comoving bolometric IR luminosity density with redshift. The right-hand axis contains the conversion 
to star-formation rate based on a scaling law by Kennicutt et al. (1998). The results from this work are presented as blue filled 
squares. The orange filled region is the results of Le Floc'h et al. (2005) up to z ~ 1, the light pink area marks the data obtained 
by Caputi et al. (2007, the z = 0.2 point comes from the bolometric IR LF derived from the 8 //m LF by Huang et al. 2007) and 
the magenta curve those by Perez-Gonzalez et al. (2005). We also report the separate contributions from LIRGs (here defined ad 
10 n L G < L]r < 1O 12 L , green dashed line) and ULIRGs (L/« > 1O 12 L , red dot-dashed line) to the IR luminosity density. The 
black dotted line corresponds to the prediction from the model of Franceschini et al. (2009). In the redshift interval < z < 1, the 
bolometric IR luminosity density evolves as (1 + Z ) 3&±0A , 



between luminosity and density evolution increases, due to the 
more limited range in luminosity covered by our LF. However, 
the LFs at z > 1 are consistent with no or limited evolution. 



7.2. The bolometric IR luminosity density 

One of the goals for computing the bolometric luminosity func- 
tions is to obtain an estimate of the total comoving IR luminos- 
ity density as a function of redshift. This inference has important 
cosmological implications, because it is tightly related with the 
evolution of the comoving Star Formation Rate density and of 
the rate of gravitational accretion into black-holes. 

Indeed, as discussed in many recent papers based mostly 
on Spitzer data (Yan et al. 2007, Sajina et al. 2007, 2008; 
Daddi et al. 2007; Fiore et al. 2008; see also Fadda et al. 2002, 
Franceschini et al. 2005), violent starburst and AGN activities 
appear often concomitant in Spitzer 24 pm sources at z ~ 2. 
All these analyses agree that, after excluding the most obvious 
type-1 AGN, it is not easy to estimate the relative AGN/starburst 
contributions in high-redshift sources, even including the widest 
SED coverage and optical/IR spectroscopic data, given the role 



of dust extinction in hiding the primary energy source and the 
very limited diagnostic power of dust re-radiated spectra. 

For this reason the only attempt we have made to differenti- 
ate the AGN contribution in our estimate of the luminosity den- 
sity was to account for the type-1 quasar contribution, whose 
power source we can entirely attribute to gravitational accretion. 
Although we might expect that a significant number of other 
among the most luminous 24 pm sources should include some 
fractional AGN contributions, it is likely that the bolometric lu- 
minosity of the bulk of the population is dominated by stellar 
processes (e.g. Ballantyne & Papovich 2007). We will in any 
case account for the AGN-dominated sources in our later analy- 
sis. 

By integrating the best-fit LF in each redshift bins, shown in 
Figure [15] we obtained our determination of the total IR lumi- 
nosity density up to z ~ 2.5, which is reported in Figure [Tolas 
blue filled squares and tabulated in Table [6] In addition to the 
Poisson noise, the errorbars that we associated to our estimates 
account for two major sources of uncertainties: the slope of the 
faint end of the luminosity function (which was set to a — 1.2 
with an error inferred from fitting the local LF), and the error on 
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the normalization parameter (<£*) in the y 1 minimization proce- 
dure. 

Within the previously mentioned uncertainties, we find a 
general agreement with previously published results based on 
IR data. The orange filled region in Figure [16] shows the re- 
sults of Le Floc'h et al. (2005) up to z ~ 1, the light pink area 
marks the data obtained by Caputi et al. (2007) and the magenta 
curve those by Perez-Gonzalez et al. (2005). In the redshift range 
< z < 1, we find that the total bolometric IR luminosity den- 
sity evolves as (1 + z) 3 8±0 - 4 . This evolution is very close to that 
derived by Le Floc'h et al. (2005), who found (1 + z) 3 - 9 between 
Z — and z = 1, but is significantly steeper than estimated by 
Caputi et al. (2007), (1 +z) 31 . At higher redshifts our results sug- 
gest a flattening of the IR luminosity density above z ~ 1 . The 
results from the IR data are also in fully agreement within what 
is derived from the far-ultraviolet data: a strong increase up to 
Z ~ 1, and a flatting above as shown for example in Tresse et al. 
(2007). 

We also report in Fig. [16] the separate contributions from 
LIRGs (here defined as objects with 10 11 < L IR < 10 12 L o, green 
dashed line) and ULIRGs (L/« > 10 12 L Q , red dot-dashed line) 
to the IR luminosity density. At z ~ 0.2, ~30% of the bolomet- 
ric IR luminosity density is contained in LIRGs and < 1% in 
ULIRGs. At z = 1, we find that LIRGs and ULIRGs contribute 
45% and 6%, respectively, to the total IR luminosity density. At 
Z ~ 2, the contributions of LIRGs and ULIRGs become 48% 
and 45% of the total budget, respectively: the cosmic evolution 
of the highest- and of the moderate-luminosity sources appear to 
be drastically different. This finding is consistent with the results 
of previous Spitzer studies of the specific SF for star-forming 
galaxies (e.g. Perez-Gonzalez et al. 2005, P. Santini et al. 2009). 

Finally, we have estimated the comoving luminosity densi- 
ties after having taken out the contribution of the 180 sources 
identified as type-1 AGNs from our SED fitting, as explained 
in Sect. [5] However, this does not make any significant differ- 
ence in Fig. [16] Assuming that the remaining part of the pop- 
ulation has a bolometric emission dominated by star-formation, 
we report on the right-hand ordinate axis of Fig.[l6]the comov- 
ing SFR density, that we calculate from Kennicutt et al. (1998): 
SFR[M Q /yr] = 1.7 x 10- 1() L([8 - I000pm])/L Q . 



If confirmed, our results are obviously much more significant 
and further reaching, because they concern the bulk of the emis- 
sion by the whole cosmic source population selected in the IR. In 
any case, this similarity of the evolutionary properties of our IR- 
selected galaxy population and those of X-ray AGNs strongly 
supports the case for co-eval AGN/starburst activity and for a 
physical relation between the two populations. 




0.0 0.5 1.0 1.5 2.0 



z 

Fig. 17. The space density of IR selected galaxies as a function 
of redshift in different luminosity classes. The numbers with the 
same colour coding at the left side of each curve indicate the 
corresponding IR luminosity in solar units (logarithmic). 



7.3. Evolution of the comoving source space density 

A more detailed inspection of the evolutionary properties of the 
IR selected sources is presented in Figure[17]in terms of the vari- 
ation with redshift of the comoving source space densities for 
different IR luminosity classes. The data points here have been 
derived from the best-fit to the bolometric LFs shown in Figure 
[T5l but we report only the data within luminosity intervals ac- 
ceptably constrained by our observations. 

The figure indicates a systematic shift with redshift of the 
number density peak as a function of luminosity: the brighest 
IR galaxies formed earlier in the cosmic history (z >~ 1.5), 
while the number denisty of the less luminous ones peaks at 
lower redshifts (z ~ 1 for L/« < 10 11 ). Although our results do 
not constrain the high-redshift evolution of the lower luminosity 
sources, they provide a tentative detection of a downsizing ef- 
fect in the evolution of the IR galaxy populations, hence in the 
cosmic star formation history. This behaviour appears to paral- 
lel a similar one reported for a completely different population 
of cosmic sources, i.e. type-1 AGNs selected in X-rays and ana- 
lyzed by Hasinger et al. (2005, see also Smolcic et al. 2009 for 
radio selected galaxies). 



7.4. Comparison with model predictions 

We report in Figures [4] [5] l9l[T2l and [16] a comparison of 
our results with predictions of a phenomenological model by 
Franceschini et al. (2009, AF09). The latter, in particular, are 
compared with the rest-frame 24 pm LF in Figure [12] (dashed 
lines), showing a generally good agreement. The low-z be- 
haviour of the AF09 model has been calibrated on a variety of 
number count data at bright fluxes (including the IRAS all-sky 
IR counts and those from the Spitzer SWIRE project), which 
guarantees excellent control of the low-redshift universal emis- 
sivity in the IR. 

Compared to previous analyses essentially based on pre- 
Spitzer surveys (e.g. Franceschini et al. 2001), an important ad- 
dition in AF09 was the introduction of an evolutionary compo- 
nent of very luminous IR galaxies dominating the cosmic IR 
emissivity at z ~2, and essentially absent or very rare locally, 
hence caracterized by an extremely fast evolution in cosmic 
time. This population is revealed from the analysis of the galaxy 
samples selected by the Spitzer/MlPS deep 24pm and SCUBA 
sub-millimetric observations, while it was undetected by ISO. 
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In Figure [16] we compare the redshift evolution of our de- 
rived cosmic star-formation rate with the AF09 model. The gen- 
eral trend is nicely reproduced, implying (together with all the 
other observational constraints considered by the model) that the 
main features of the IR galaxy evolution are understood. Our re- 
sults support a scenario where the fast evolution observed from 
z — to 1 flattens above z ~ 1 and keeps approximately flat up 
to z ~ 2.5. Moreover, both models and the present data provide 
clear evidence for the existence of a population of very luminous 
galaxies becoming increasingly important at z > 1 . The current 
interpretation identifies these objects with the progenitors of the 
spheroidal galaxies. A more detailed discussion is deferred to 
AF09. 

8. Summary 

In this paper we have exploited a combination of data from the 
GOODS and VVDS-SWIRE multi-wavelength surveys to deter- 
mine mid-IR and bolometric luminosity functions (and thus to 
estimate the SFR density) over a wide (0 < z < 2.5) redshift 
interval. 

The primary selection for this analysis comes from flux- 
limited samples of MTPS/Spitzer 24 fim sources in the GOODS 
(North and South) and VVDS-SWIRE fields, at the flux lim- 
its of 5(24/im) >80 /zJy and S(24fim) >400 //Jy, respectively. 
Our combined use of these two sensitivity thresholds helped us 
to obtain a wide coverage of the luminosity-redshift plane with 
good statistics. We have performed a careful identification of the 
optical counterparts, by applying a maximum likelihood tech- 
nique for the deeper GOODS surveys, while the identification 
was more straightforward at the brighter VVDS-SWIRE limits. 
We took advantage of the extensive spectroscopic information 
available in the GOODS fields and of the high-quality photo- 
metric redshifts in the VVDS-SWIRE area, combined with the 
rich multiwavelength optical-to-IR datasets, to perform a fitting 
procedure on the whole SED for each source in our sample. In 
this way, we could estimate the rest-frame luminosity at differ- 
ent wavelengths, at the same time obtaining the spectral extrap- 
olations needed to compute the bolometric luminosity for each 
object. 

We have then computed rest-frame luminosity functions at 
8, 12, 15 and 24 /vm for comparison with a variety of previously 
published results at these wavelengths. We found a fairly good 
agreement with previous results at the various rest-frame wave- 
lengths and redshift bins. 

We extrapolated total IR ([8-1000 fim]) luminosities from 
our best-fit to the SEDs of each source, and used these to derive 
the bolometric LF up to z ~ 2.5. Adopting our fitting function, 
the number density of sources, as parameterized by O*, evolves 
as e (1 + z) 1 1 in the redshift range < z < 1, while the typi- 
cal infrared luminosity (L*) evolves as (1 + z) 2J '. Above z > 1, 
the evolution degeneracy between number density and luminos- 
ity is more critical, but the comoving density is bound to keep 
roughly constant with redshift. By integrating the best-fit LF in 
each redshift bin, we obtained a robust determination of the IR 
luminosity density up to z ~ 2 and we found a general agree- 
ment with previously published results based on IR data. In the 
redshift range < z < 1, we find that the bolometric IR luminos- 
ity density evolves as (1 + z) 38±04 . At higher redshifts our result 
seem to confirm a flattening of the IR luminosity density. 

We estimated the separate contributions from LIRGs (10 11 < 
L IR < 1O 12 L ) and ULIRGs (L IR > 1O 12 L ) to the IR luminos- 
ity density. At z ~ 0.2, ~30% of the bolometric IR luminosity 
density comes from LIRGs and < 1% from ULIRGs. At z = 1, 



the fractional contributions from LIRGs and ULIRGs are 45% 
and 6%, respectively. At z ~ 2, the contributions of LIRGs and 
ULIRGs become 48% and 45% of the total budget, respectively. 
This demonstrates that the cosmic evolution of the highest- and 
of the moderate-luminosity sources are therefore drastically dif- 
ferent. 

Also in comparison with model predictions by Franceschini 
et al. (2009, in preparation, but see also Franceschini et al. 2008), 
our results confirm the presence of a very rapid increase of the 
galaxy IR volume emissivity up to z ~ 1 and the evidence for the 
existence of a population of very luminous galaxies becoming 
dominant at z > 1 . 

Finally, our data seem to indicate a clear dependence on lu- 
minosity of the source comoving number density peak as a func- 
tion of redshift: the brighest IR galaxies formed stars earlier in 
the cosmic history (z >~ 1.5), while the star-formation activ- 
ity for the less luminous arises later (z ~ 1 for Ljr < 10 11 ). 
This confirms a downsizing pattern in the evolution of the IR 
emissivity of galaxy populations. This behaviour is at least qual- 
itatively similar to what has been already detected for AGNs at 
other wavelengths, although our limited coverage of the faint lu- 
minosity end prevents us from achieving more definite conclu- 
sions. 
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Table 1. 8 /J.m luminosity function values. The units are log(O) [Mpc 3 dex ']. 
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Table 2. 12 //m luminosity function values. The units are log(cD) [Mpc 3 dex ']. 
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Table 3. 15 //m luminosity function values. The units are log(<l>) [Mpc 3 dex ']. 
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Table 4. 24 //m luminosity function values. The units are log(<l>) [Mpc 3 dex ']. 
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Table 5. Total IR (bolometric) luminosity function values. The units are log(cD) [Mpc 3 dex ']. 
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Table 6. Bolometric IR luminosity density values as a function 
of redshift. 



